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Abstract. 

Solutions to the field equations generated from Lagrangians of the form f(R) are 
considered. The spherically symmetric solutions to these equations are discussed, 
paying particular attention to features that differ from the standard Schwarzschild 
solution. The asymptotic form of solutions will be described, as will the lack of validity 
of Birkhoff 's theorem. Exact solutions are presented which illustrate these points and 
the stability and equations of motion of particles in these space-times are investigated. 
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1. Introduction 

It has been known for some time that Birkhoff 's theorem and asymptotic flatness are not 
ensured in the vacuum solutions of generalised fourth-order theories of gravity, as they 
are in general relativity (GR) (see e.g. [HE])- The standard approach to overcoming 
this problem is to consider theories which contain in their generating Lagrangian an 
Einstein-Hilbert term which dominates the field equations in the low-curvature limit. 
In this way the usual behaviour of approaching Minkowski space at asymptotically 
large distances from sources is assumed to occur and perturbative expansions about a 
Minkowski background are then possible. The weak-field limit of the theory in question 
can then be investigated in a straightforward way and compared with experiment. 

Whilst the existence of a Minkowski background is a great simplification, and 
very useful in investigations of the weak-field, it is unclear whether or not a theory 
should be disregarded solely on the grounds of this limit not existing. There are a 
number of gravitational theories that can be conceived of which do not always admit 
Minkowski space as a solution. These include theories derived from Lagrangians of the 
form R 1+s (where 5^0)[3lHIHlEllZllH!,as well as those derived form R + a/R n 
(where n > 0) pm EH] ■ The latter has caused considerable debate as to whether or not 
it is compatible with solar system experiments. Studies on this subject usually follow 
the prescribed analysis for computing the weak field geometry in GR; picking a highly 
symmetric background^] and calculating the form of spherically symmetric perturbations 
by linearising the field equations and solving them, to first order in perturbations. 
However, it is often unclear how an appropriate background to expand about should be 
chosen, and how the lack of Birkhoff 's theorem or asymptotic flatness should be taken 
into account when performing these analyses. We will discuss these points here and 
attempt to make some progress into understanding the spherically symmetric situation 
by using explicit exact vacuum solutions of the field equations. 

In section 2 we present the field equations for the theories we will be considering 
and two exact spherically symmetric solutions to these field equations. In section 3 we 
perform a linear perturbation analysis about the backgrounds of the two exact solutions 
and in section 4 we find the equations of motion for test-particles in these space-times. 
In section 5 we summarise our results. 

2. Field equations and exact solutions 

We will consider here gravitational theories derived from a Lagrangian of the form 

C = f(R) (1) 

where f(R) is an arbitrary power-series expandable function of the Ricci curvature scalar 
R. Extremizing the action that is obtained by integrating ([1]) over all space then gives 

| often Minkowski space, even though it is not always a solution of the field equations 
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the field equations 



f'Rab ~ ^f'9ab + f'- Cd {9ab9cd ~ 9ac9bd) = 8nGT ab (2) 

where matter sources have been included, primes denote differentiation with respect to R 
and we have disregarded boundary terms. It can be seen that in the low-curvature limit, 
R << 1 (assuming this limit exists), the term in the power-series expansion of f(R) with 
the lowest power will be the one which dominates the field equations. For this reason 
we specialise our considerations, from this point on, to the choice f(R) = R 1+s which 
reduces to general relativity in the limit 5 — > 0. This choice is the theory considered in 
[HI SI [5] and will be the effective theory at low curvatures for f(R) = R + a/R n when 
6 = — (1 + n). The field equations are then given by 

5(1 - 5 2 )R S ^± - 5(1 + 5)R S ^ + (1 + 6)R 5 R ab - ±g ab RR s (3) 

- g ab 6(l - 5 2 )R 5 ^- + 5(1 + 5)g ab R 5 ™ = 8vrT ab . 

We will now present two exact spherically symmetric vacuum solutions to this set of 
equations. The first of these was found previously in [3] and is given by the line-element 

dv 

Solution 1 ds 2 = -A l (r)dt 2 + — — + r W (4) 



where 



Bi(r) 



25 ^1 , Cl 

5l(r) = (1 -25 + A5 2 )(l-25(l + 5)) ' ' 



(l-2^+4(S 2 ) 

r 



and C% is a constant. This solution is conformally related to the Q = limit of a solution 
found by Chan, Home and Mann [12] and reduces to the Schwarzschild solution in the 
limit 5^0. The second solution is given by 

Solution 2 ds 2 = -A 2 (r)dt 2 + a 2 (t)B 2 (r)(dr 2 + r 2 dVl 2 ) (5) 

where 



Mr)= — & a(t)=t°J 




s (1+26) 



B 2 (r)= 1 + — A^ 26 - 1 q 2 = 1-25 + 45 



and C*2 is a constant. This solution is conformally related to one found by Fonarev 
and again reduces to the Schwarzschild solution in the limit 5^0. 

These two solutions can be seen to exhibit features not present in the Schwarzschild 
solution of GR. Both of these solutions are strongly curved, but each displays this 
curvature in a different way. Solution 1 is static and does not reduce to an r independent 
form in the limit r — > oo (despite the Ricci scalar approaching in this limit). Solution 
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2 displays more conventional behaviour in the limit r — > oo, but shows explicitly the 
lack of validity of Birkhoff 's theorem. This solution becomes r independent in the limit 
r — > oo, but still displays strong curvature in this limit as the metric reduces to the 
spatially flat vacuum Friedmann- Robert son- Walker cosmology found in [3]. 

We will now continue to find the general form of spherically symmetric 
perturbations to the backgrounds (r — > oo limit) of the two exact solutions above. 
It will be seen that there exist extra modes which are not excited in the exact solutions, 
but that the modes corresponding to the linearised exact solutions above are the ones 
most interesting for performing gravitational experiments in these space-times. 



3. Linear perturbation analysis 

Perturbative analyses in the literature are often performed about Minkowski space or 
de-Sitter space. This is perfectly acceptable practise in GR and fourth-order theories in 
which an Einstein-Hilbert term dominates in the low curvature regime. In other fourth- 
order theories, of the type considered here, in which the Einstein-Hilbert term does not 
dominate the low curvature regime then there is good reason to consider perturbing 
about other backgrounds. We have shown explicitly, with exact solutions, the existence 
of other spherically symmetric space-times. We will now proceed to perform a linear 
perturbation analysis about the backgrounds (r — > oo limit) of these two space-times. 
The general solution to first order in perturbations will be found and it will be shown 
that the form of these linearised solutions will be strongly dependant on the background. 



3.1. Solution 1 

An analysis of the linear perturbations about the background prescribed by equation 
P| has already been performed in [3]. We will simply quote the result of this study 
here, the reader is referred to [3] for details of the derivation. 
Writing the perturbed line-element as 

ds 2 = -r 25 l^T(l + V(r))dt 2 

(1 -2<5 + 4<5 2 )(l-25-2<5 2 ) M xxr , ... 2 2jo2 

+ 09 -(1 + W(r))dr 2 + r 2 dtt 2 , 

(1 -b) 1 

substituting this into the field equations and linearising them to first order in V and W 
allows the general solution to first order in perturbations to be found as 

V(r) = ciV\(r) + c 2 V 2 (r) + c 3 V 3 (r) + constant 
W{r) = -ciVi(r) + c 2 W 2 (r) + c 3 W 3 {r) 



where 



(1-25+4S 2 ) 

Vi = — r 

V 2 = 2(2 l _ m2 - m3) ((1 + 26) 2 sin(Alogr) + 2A(1 - 6) cos(Alogr)) 
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W 2 = r 2(i=«i sin(Alogr) 

(1-2S+4S 2 ) 



V? 



1 + 25) ((1 + 25) 2 cos(Alogr) - 2A(1 - 5) sin(A logr)) r 



2(1-4) 
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2(2 - 35 + 125 2 + 165 3 ) 



and 



(l-2<5+4<5 2 ) 

W 3 = r 2 (!- a ) cos(Alogr) 



V7 - 285 + 365 2 - 165 3 - 805 4 
A = 



2(1-5) 

The extra constant in V can be absorbed into the definition of the time coordinate. This 
procedure of solving the linearised field equations does not ensure that the solution 
obtained will be the linearisation of the general solution to the full non-linear field 
equations, but it is encouraging to note that the c 2 mode corresponds to the linearisation 
of the exact solution flU). 



3.2. Solution 2 

The corresponding perturbative analysis about the background given by (jSJ) will now be 
performed. Writing the line-element as 

ds 2 = -(1 + P(r))dt 2 + b 2 (t)(l + Q(r))(dr 2 + r 2 dVl 2 ) (6) 

allows the vacuum field equations to be linearised in P and Q. These linearised field 
equations are given in the appendix, and are solved to give the solutions 

p = _e, 2c 5 (l- 6S + 4P +4*3) 

r (5 - 145 - 125 2 ) 

^ (l-25)c 4 . 2 

Q = h oc 5 r + constant 

r 

where C4 and C5 are constants and the two other constant terms in P and Q are 
independent of each other and can be absorbed into t and s by redefinitions. Again, it 
can be seen that one of the modes, C4, corresponds to the linearised version of the exact 
solution, (jSJ). 

It can immediately be seen that the form of these linearised solutions are quite 
different to those obtained by expanding around the background of solution 1, (j3J). 
Whist the expansion about (j4j) produces damped oscillatory modes, as well as the mode 
corresponding to the linearised exact solution, the expansion about (jSJ) produces more 
familiar looking terms proportional to r 2 . Aside from the different form of these extra 
modes there is also a noticeable difference in the terms corresponding to the linearised 
exact solutions, which both go as r _1 in the limit 5^0, but behave differently from each 
other when 5^0. This shows explicitly the differences that can arise when linearising 
about different backgrounds. Not only are there extra modes which can take different 
functional forms, but even the modes which reduce to the Schwarzschild limit as 5 — > 
are different, depending on the background that has been chosen. 



Spherically Symmetric Solutions to Fourth- Order Theories of Gravity 



6 



We will now proceed to calculate the equations of motion of particles following 
geodesies of these space-times, to Post-Newtonian order. It will be shown that not 
only are the terms due to the linear perturbations different (as should be expected as 
the perturbations themselves have been shown to be background dependant) but that 
the background itself also contributes extra terms to the post-Newtonian equations of 
motion. 



4. Equations of motion 

We will now proceed to calculate the geodesies of the two space-times, (jlj) and In 
doing this we will neglect the contributions from the oscillating modes to (jlj) and the 
contribution of the r 2 mode to (jSJ) so that we only take into account the modes which 
go as r _1 , in the limit 5 — > 0. These are the modes corresponding to the linearisation of 
the exact solutions (jlj) and (jol) . In performing this computation we will transform the 
solution (jlj) into isotropic coordinates (details of this transformation are given in [3]). 
The geodesic equation can be written, as usual, in the form 

L = 

dX 2 13 dX dX 

where A can be taken as proper time for a time-like geodesic, or as an affine parameter 
for a null geodesic. In terms of coordinate time this can be re-written as 

d'.rH + f _ r o ^\ chfcW = Q 



dt 2 \ 13 13 dt J dt dt 
Substituting the linearised solutions into this equation will then give the equations of 
motion for test particles in these space-times. 



4-1. Solution 1 

In isotropic coordinates the linearised version of (jlj) can be written as 
ds 2 = -A 3 (r)dt 2 + B 3 (r)(dr 2 + r 2 dQ 2 ) 



where 



2£(l ± 2 i ) / (1 - 2<5)c 6 

A 3 (r) — r « II 



r 

B 3 (r) = r n 1 -| 



ni 



and 



r 



n 2 = (1 - 28 - 28 Z ){1 - 25 + AS' 
m 2 



[1-25 + A5 2 ) 



(1 - 25 - 25 2 ) 

and c 6 is constant. Substituting this metric into the geodesic equation (j7|) gives, to 
post-Newtonian order, 

d 2 x 5(1 + 25) (1 - 85 + A5 2 )c 6 (l-8^ + 4^ 2 )c 2 

~~JpI ~ A_i er ~ 3(1-2<S) e r H 4(1-^)2 1 ^ 

at nrm 2nr — s 1 2nr^H — 1 
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[1 — 5 — n) mcQ \ / dx 
r 1 ~ 2r 1+m / V ~dt. °' 



+ 



2(m — 1) 2(1— 5)mcQ\ dxdx 



1+m / dt dt 



which to first order in S is 
d 2 x Gm 



dt 2 r 2 



dt 



* G 2 m 2 Gm dxdx 5 

e r + 4 — — e r + A—-e r ■ — — - -e r (9) 

r d at at r 



where the Newtonian limit has been used to set Cg. It can be seen that these equations 
of motion are modified from their usual form in GR, both in the pre-multiplicative 
factors of the terms with GR counterparts as well as modifications to the powers of r 
and the existence of entirely new terms which vanish in the limit 5^0. The first order 
corrections due to small, but non-zero, 5 are primarily due to the new term arising in 
([§]). This term corresponds to a new force which drops off as r _1 and is the term which 
was used in [3] to impose upon 5 the tight constraint 

5 = 2.7 ±4.5 x KT 19 (10) 

from observations of the perihelion precession of Mercury [TJ]. 

4-2. Solution 2 

Solution 2 has already been given in isotropic coordinates in equation ([51), which on 
substitution into ((Tj) gives to post-Newtonian order the equation of motion 

g? 2 x Gm , r .G 2 m 2 nr .Gm ( dx\ 2 ,„ . 

e r + 4(l-5)^— e r -(l-25)— 5- — (11) 



where 



_ N Gm dx dx „ , dx 

+ 4(1 - 5)— ^-e r • — H— 

r z dr dr dr 



Ida 
At = - — , 

a dr 




the conformal time coordinate r is defined by dt = adr and C4 has been set by the 
appropriate Newtonian limit. The equations of motion for this solution are considerably 
simpler than those of solution 1, but still differ from those of GR in significant ways. All 
of the terms except the last in equation ( fill have GR counterparts and the powers of r 
in these terms are all the same as in the general relativistic case. The premultiplicative 
factors of these terms are, however, modified and can be described adequately within 
the frame-work of the parameterised post-Newtonian approach [15] by assigning 

13=1 and 7 = 1 - 26. 

By making this identification the constraints on 7 from observations of the Shapiro time 
delay of radio signals from the Cassini space probe [IE] can be used to impose upon 5 
the constraint 

5 = -1.1 ±1.2 x KT 5 . (12) 
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As well as the usual effects associated with 7 — 1 being non-zero there are extra effects 
in this space-time due to the last term in pip . This term is proportional to the velocity 
of the test-particle (when v « c) and is zero for photons. For this reason we identify 
it as a friction term, with the friction coefficient being given by 7i. This 'friction' is 
a purely gravitational effect and is not due to any non-gravitational interaction of test 
particles with any other matter. 

5. Conclusions 

We have considered here the problem of finding spherically symmetric vacuum solutions 
of fourth-order theories of gravity. In these theories the Einstein-Hilbert term may not 
be the dominating one at low curvatures in the generating Lagrangian and the validity 
of taking a Minkowski or de-Sitter background to expand about cannot be assumed so 
readily. As well as this, Birkhoff's theorem is not valid in fourth-order theories, which 
allows for a much larger manifold of solutions, when spherical symmetry is imposed. 
These features allow for the interesting possibility of having a choice of backgrounds 
to expand about, both static and non-static, which can vary considerably in form from 
their general relativistic counterparts. 

We have presented two exact solutions which are spherically symmetric and reduce 
to the Schwarzschild solution in the limit of the fourth-order theory reducing to general 
relativity. These solutions have very different backgrounds (r — > 00 limits) and 
gravitational physics in each of them is correspondingly different. The first of these 
solutions is static and has a non-trivial dependence on the radial coordinate r, as r — > 00. 
The second is non-static and has an r independent form in the limit r — > 00. In some 
sense these solutions can be considered as being extremes of the general solution; in the 
first case the strong curvature at large distances from the centre of symmetry is only 
dependent on r and in the second the curvature at large distances is only a function of 
t (in the limit r — > 00). 

Using these exact solutions we have shown how different choices of background lead 
to very different solutions of the linearised field equations, when spherically symmetric 
perturbations are introduced. By taking the perturbation modes which reduce to 
something approximating Newtonian gravity, in the appropriate limit, we have derived 
the equations of motion for test particles following geodesies of these space times. The 
equations of motion are found to be strongly dependent on the choice of background. 
In both cases the form of the terms that have counterparts in GR are dependent on 
the background. As well as this dependence there are additional terms due to the 
background itself which have no counterparts in the general relativistic solutions. For 
the static case these extra terms look like a fifth force which drops off as ~ r _1 and for 
the time dependent case they look like a friction term oc v, when v << c. 

These considerations show explicitly that there are a number of different 
backgrounds about which one may choose to perform a perturbative expansion, in 
fourth-order theories, and that these backgrounds can display behaviour which is 
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not permitted in general relativity. Moreover, it has been shown that the choice of 
background is highly non-trivial when performing a perturbative expansion. The form of 
the perturbations are background dependent and, correspondingly, so are the geodesies 
of the perturbed space-times. 

The existence of multiple backgrounds, all of which appear to be valid solutions of 
the field equations, prompts the question: which is the appropriate solution for physical 
situations? We will not attempt to address this question fully here, but will offer a 
few suggestions and leave a more comprehensive analysis for a future study. Firstly, 
a physically relevant background solution should be stable to perturbations. Imposing 
the extra symmetry of time-independence it was shown in [3] that the r — > oo limit of 
solution (j3J) is the attractor of the general spherically symmetric and time independent 
solution, in this limit. It was also shown in reference [3] that the r — > oo limit of solution 
(J5j) is the attractor of the general homogeneous and isotropic vacuum solution in the 
limit t — > oo (at least for some range of 5). In this sense both of these solutions can be 
considered as stable, but this does not allow us to differentiate between the two of them 
which is the physically more relevant. In order to establish the appropriate solution 
for any realistic physical situation it will be necessary to apply boundary conditions to 
either or both of an interior solution for an energy-momentum distribution and/or the 
relevant cosmological solution. Whilst a number of cosmological solutions are known 
[31 HI El E] this problem is made more difficult by a lack of exact solutions which may 
be valid in the interior of, for example, a star. A more comprehensive analysis is also 
required in order to establish the effect of allowing the perturbations themselves to be 
time dependent. 

6. Appendix 

Substituting the perturbed line-element (jHJ) into the vacuum field equations ([3]) allows 
the following set of coupled, linearised differential equations to be found. 
The t — t equation is 



65(1 - 25 - 25 2 )(5 - 145 - 125 2 )VP 

-12(1 - 25 - 25 2 )(1 - 65 + 45 2 + 45 3 )VQ 

_ 2(1-2S-2S 2 ) 

+(1 + 25) (1 - 5) 2 t (i-*) = 



the r — r equation is 

(2 - 175 + 185 2 + 525 3 + 85 4 )P" -2(1-5 + 35 2 



345 3 - 325 4 ) — 



r 



+2(2 - 85 - 95 2 + 165 3 + 85 4 )Q" + 2(1 - 45 



155 2 + 205 3 + 165 4 ) 



r 




and the t — r equation is 

ip' = 
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where 

■0 = VP + 2VQ, 

primes here denote differentiation with respect to r and V is the Laplacian on the three 
dimensional subspace. The 9 — 9 and — equations are then linear combinations of 
the above equations. 
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